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ABSTRACT 

Based on a new approach on modeling the magnetically dominated outflows from AGNs (Li et al. 2006), we 
study the propagation of magnetic tower jets in gravitationally stratified atmospheres (such as a galaxy cluster 
environment) in large scales (> tens of kpc) by performing three-dimensional magnetohydrodynamic (MHD) 
simulations. We present the detailed analysis of the MHD waves, the cylindrical radial force balance, and the 
collimation of magnetic tower jets. As magnetic energy is injected into a small central volume over a finite 
amount of time, the magnetic fields expand down the background density gradient, forming a collimated jet and 
an expanded "lobe" due to the gradually decreasing background density and pressure. Both the jet and lobes 
are magnetically dominated. In addition, the injection and expansion produce a hydrodynamic shock wave 
that is moving ahead of and enclosing the magnetic tower jet. This shock can eventually break the hydrostatic 
equilibrium in the ambient medium and cause a global gravitational contraction. This contraction produces a 
strong compression at the head of the magnetic tower front and helps to collimate radially to produce a slender- 
shaped jet. At the outer edge of the jet, the magnetic pressure is balanced by the background (modified) gas 
pressure, without any significant contribution from the hoop stress. On the other hand, along the central axis of 
the jet, hoop stress is the dominant force in shaping the central collimation of the poloidal current. The system, 
which possesses a highly wound helical magnetic configuration, never quite reaches a force-free equilibrium 
state though the evolution becomes much slower at late stages. The simulations were performed without any 
initial perturbations so the overall structures of the jet remain mostly axisymmetric. 

Subject headings: magnetic fields — galaxies: active — galaxies: jets — methods: numerical — magnetohy- 
drodynamics (MHD) 



1. INTRODUCTION 

A number of astronomical systems have been discovered 
to eject tightly collimated and hyper-sonic plasma beams and 
large amounts of magnetic fields into the interstellar, intr- 
acluster and intergalactic medium from the central objects 
during their initial/final (often violent) stages. Magnetohy- 
drodynamic (MHD) mechanisms are frequently invoked to 
model the launching, acceleration and collimation of jets 
from Young Stellar Objects (YSOs), X-ray binaries (XRBs), 
Active Galactic Nuclei (AGN s), Microquasars, an d Quasars 
(QSOs) (see, e. g. . iMeier. Koide. & Uchidal feOOl. and refer- 
ences therein). 

Theory of magnetically driven out flows in the electro mag- 
netic regime ha s been proposed by Blandford ( 1976|) and 
lLovelacel \ 19761) and subsequentl y applied to rotating black 
holes ( Blandford & Znai ekll 19771) and to magnetized accre- 
tion disks ( Blandford & Pavne 1982). By definition, these 
outflows initially are dominated by electromagnetic forces 
close to the central engine. In these and subsequent models of 
magnetically driven outflows (jets/winds), the plasma veloc- 
ity passes successively through the hydrodynamic (HD) sonic, 
slow-magnetosonic, Alfvenic, and fast-magnetosonic critical 
surfaces. 

The first attempt to investigate the nonlinear (time- 
dependent) behavior of magnetic ally driven outflows from 
accretion disks was performed bv lUchida & Shibatal d!985h . 
The differential rotation in the system (central star/black hole 
and the accretion disk) creates a magnetic coil that simulta- 
neously expels and pinches some of the infalling material. 
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The buildup of the azimuthal (toroidal) field component in 
the accretion disk is released along the poloidal field lines as 
large-amplitude torsional Alfven waves ("sweeping magnetic 
twist"). After their pioneering work, a number of numerical 
simulations to study the MHD jets have been done (see, e.g., 
Ferrari! [19981 and references therein). An underlying large- 
scale poloidal field for producing the magnetically driven jets 
is almost universally assumed in many theoretical/numerical 
models. However, the origin and existence of such a galactic 
magnetic field are still poorly understood. 

In contrast with the l arge-sc ale field models , Lynd en-Bell 
llLvnden-Beil&Boilvl fl99l iLvnden-Belll 119961 12001 
2006) examined the expansion of the local force-free mag- 
netic loops anchored to the star and the accretion disk 
by using the semi-analytic approach. Twisted magnetic 
fluxes due to the disk rotation make the magnetic loops 
unstable and splay out at a semi-angle 60° from the rota- 
tional axis of the disk. Global magnetostatic solutions of 
magnetic towe rs with externa l thermal pressure were also 
computed by iLi et al. I J2001I) usin g the Grad-Shafranoy 
equation in axisymmet ry (see al so. lLov elace et akl 12002): 
Lovelace & Romanova 2003; Uzdenskv & MacFadven 
2006). Full MHD numerical simulations of magnetic towers 
have been performed in two-dimension (axisymmetric) 

| 1999t lUstvugova et al. 1120001 iKudoh M atsumoto & Shibata 
I2OO2I iKato. Havashi. & Matsumotol 12004 and three- 
dimension (Kato. Mineshige. & Shibata 2004). Magnetic 
towers are also observed in the labora tory experiments 
dHsu & Bellan l2002tlLebedev et al. 120051) . 

This paper describes the nonlinear dynamics of propagat- 
ing magnetic tower jets in large scales (> tens of kpc) based 
on three-dimensional MHD simulations. We follow closely 
the approach described in Li et al. (2006; hereafter Paper I). 
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Different from Paper I, which studied the dynamics of mag- 
netic field evolution in a uniform background medium, we 
present results on the injection and the subsequent expansion 
of magnetic fields in a stably stratified background medium 
that is described by an iso-thermal King model llKing 119621) . 
Since the simulated magnetic structures traverse several scale 
heights of the background medium, we regard that our simula- 
tions can be compared with the radio sources inside the galaxy 
cluster core regions. Due to limited numerical dynamic range, 
however, the injection region (see Paper I) assumed in this 
paper will be large (a few kpc). Our goal here is to provide 
the detailed analysis of the magnetic tower jets, in terms of 
its MHD wave structures, its cylindrical radial force balance, 
and collimation. The paper is organized as follows: In §2, we 
outline the model and numerical methods. In §3, we describe 
the simulation results. Discussions and conclusions are given 
in §4 and §5. 

2. MODEL ASSUMPTIONS AND NUMERICAL METHODS 

The basic model assumptions and numerical treatments we 
adopt here are essentially the same as those in Paper I. Mag- 
netic fluxes and energy are injected into a characteristic cen- 
tral volume over a finite duration. The injected fluxes are not 
force-free so that Lorentz forces cause them to expand, in- 
teracting with the background medium. For the sake of com- 
pleteness, we show the basic equations and other essential nu- 
merical setup here again, and refer readers to Paper I for more 
details. 

2.1. Basic Equations 

We solve the nonlinear system of time-dependent ideal 
MHD equations numerically in a 3-D Cartesian coordinate 
system (x, y, z): 



^+V-(pV)=p i n j (1) 

V- (pW +p + B 2 /2-BB)=-pVip (2) 
dE 

— + V • [(E + p + B 2 /2) v B(v B)] = - P V ■ S7tf> 
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V x (V x B) = B i; 



(3) 
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Here p, p, V, B, and E denote the mass density, hydrody- 
namic (gas) pressure, fluid velocity, magnetic field, and to- 
tal energy, respectively. The total energy E is defined as 
E = p/(l - 1) + pV 2 /2 + B 2 /2, where 7 is the ratio of 
the specific heats (a value of 5/3 is used). The Newtonian 
gravity is — V^. Quantities pi n j, B^, and B; n j represent the 
time-dependent injections of mass, magnetic flux, and energy, 
whose expressions are given in Paper I. 

We normalize physical quantities with the unit length scale 
i?o, density po, velocity Vb in the system, and other quanti- 
ties derived from their combinations, e.g., time as RqIVq, etc. 
These normalizing factors are summarized in Tabled Here- 
after, we will use the normalized variables throughout the pa- 
per. Note that a factor of 4ir has been absorbed into the scal- 
ing for both the magnetic field B and the corresponding cur- 
rent density J. To put our normalized physical quantities in 
an astrophysical context, we use the parameters derived from 
the X-ray observations of the Perseus cluster as an example 
(Churazov et al. 2003). These values are also given in Table 
□ 



The system of dimensionle ss equations is in tegrated in time 
by using an upwind scheme tLi & Li| |2003). Computations 
were performed on the parallel Linux clusters at LANL. 

2.2. Initial and Boundary Conditions 

One key difference from Paper I is that we now introduce 
a non-uniform background medium. An iso-thermal model 
iKing I n~962h has been adopted to model a gravitationally 
stratified ambient medium. This is applicable, for example, 
to modeling the magnetic towers from AGNs in galaxy clus- 
ters. 

The initial distributions of the background density p and gas 
pressure p are assumed to be 



p = p = 



1 



(5) 



y 2 + z 2 ) 1 ! 2 is the spherical radius and 



where R = (x 2 
R c the cluster core radius. (In the following discussion, both 
"transverse" and "radial" have the same meaning, referring 
to the cylindrical radial direction.) The parameter k controls 
the gradient of the ambient medium. Furthermore, we assume 
that the ambient gas is initially in a hydrostatic equilibrium 
under a fixed (in time and space) yet distributed gravitational 
field —Vip(R) (such as that generated by a dark matter poten- 
tial). From the initial equilibrium, we get 
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In the present paper, we choose R c — 4.0 and k = 1.0. The 
initial sound speed in the system is constant, C s |t = o = 7 1 ' 2 ~ 
1.29, throughout the computational domain. An important 
time scale is the sound crossing time r(= R/C$) w 0.78, 
normalized with r s o(= Rq/C s q) w 10.0 Myrs. Therefore, a 
unit time scale t = 1 corresponds to 12.8 Myrs. 

The total computational domain is taken to be \x\ < 16, 
\y\ < 16, and \z\ < 16 corresponding to a (80 kpc) 3 box 
in the actual length scales. The numbers of grid points in 
the simulations reported here are N x x N y x N z — 240 x 
240 x 240, where the grid points are assigned uniformly in 
the x, y, and z directions. A cell Ax (= Ay — Az ~ 0.13) 
corresponds to ~ 0.65 kpc. We use the outflow boundary 
conditions at all outer boundaries. Note that for most of the 
simulation duration, the waves and magnetic fields stay within 
the simulation box, and all magnetic fields are self-sustained 
by their internal currents. 

2.3. Injections of Magnetic Flux, Mass, and Energy 

The injections of magnetic flux, mass and its associated en- 
ergies are the same as those described in Paper I. The ratio 
between the toroidal to poloidal fluxes of the injected fields is 
characterized by a parameter a — 15, which corresponds to 
~ 6 times more toroidal flux than poloidal flux. The magnetic 
field injection rate is described by 7& and is set to be 7& = 3. 
The mass is injected at a rate of j p = 0.1 over a central vol- 
ume with a characteristic radius r p = 0.5. Magnetic fields 
and mass are continuously injected for t- m j — 3.1, after which 
the injection is turned off. These parameters correspond to 
an magnetic energy injection rate of ~ 10 43 ergs/s, a mass 
injection rate of ~ O.O46Af /yr, and an injection time ~ 40 
Myrs. 

In summary, we have set up an initial stratified cluster 
medium which is in a hydrostatic equilibrium. The magnetic 
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flux and the mass are steadily injected in a central small vol- 
ume with a radius of 1 . Since these magnetic fields are not in 
a force-free equilibrium, they will evolve, forming a magnetic 
tower and interacting with the ambient medium. 

3. RESULTS 

In this section we examine the nonlinear evolution and the 
properties of magnetic tower jets in the gravitationally strati- 
fied atmosphere. 

3.1. Overview of Formation and Propagation of a Magnetic 

Tower Jet 

Before considering our numerical results in detail, it is in- 
structive to give a brief overview of the time development of 
the magnetic tower jet system. We achieve this by present- 
ing the selected physical quantities using two-dimensional 
x — z slices at y = 0. The distributions of density at var- 
ious times (t = 2.5, 5.0, 7.5, and 10.0) are shown in Fig. 
G] At the final stage (t = 10.0), we see the formation of a 
quasi-axisymmetric (around the jet axis) magnetic tower jet 
with low density cavities (a factor of ~ 30 smaller than the 
peak density). Inside these cavities, the Alfven speed is large 
V A > 5.0, while plasma /3 is small (/3 = 2p/B 2 < 0.1). 
However, the temperature T (oc ) becomes large T > 2.5 
(~ twice the initial constant value); that is, the hotter gas is 
confined in the low-/? magnetic tower. The jet possesses a 
slender hourglass-shaped structure with a radially confined 
"body" for \z\ < 3, which is likely due to the background 
pressure profile having a core radius R c — 4. As the magnetic 
tower moves into an increasingly lower pressure background, 
the expanded "lobes" are formed. A quasi-spherical hydrody- 
namic (HD) shock wave front moves ahead of the magnetic 
tower. 

A snapshot of the gas pressure change ratio Ap/pi = 
p/pi — 1 (where p\ represents p\t=o) at t = 10.0 is shown 
in Fig. [2] Positive Ap can be seen at both the post-shock re- 
gion of the propagating HD shock wave and just ahead of the 
magnetic tower Qz\ ~ 8 — 10). The distribution of Ap forms 
a "U"-shaped bow-shock-like structure around the head of the 
magnetic tower. This structure, however, does not appear un- 
til t ~ 7.5. This is apparently caused by the local compres- 
sion between the head of the magnetic tower and the reverse 
MHD slow-mode wave (see discussions in the next section). 
The gas pressure inside the magnetic tower becomes small 
(\Ap/pi\ < 0.5) due to the magnetic flux expansion. Note 
that the light-blue region between the HD shock and the mag- 
netic tower shows a small pressure decrease (Ap/pi « —0.1). 
In later sections, we will discuss the origin of this depression 
and what role it plays in the dynamics of magnetic tower jets. 

The magnetic tower jet has a well-ordered helical magnetic 
configuration. The 3-D view of the selected magnetic lines of 
force, as illustrated in Fig. [3] indicates that a tightly wound 
central helix goes up along the central axis and a loosely 
wound helix comes back at the outer edge of the magnetic 
tower jet. The magnetic pitch B^ / y/ B% + B£ has a broad 
distribution with a maximum of ~ 15. Figure|4]shows a snap- 
shot of the axial current density J z at t = 10.0. Clearly, the 
axial current flow displays a closed circulating current system 
in which it flows along the central axis (the "forward" current 
J F ) and returns on the conically shaped path that is on the 
outside (the "return" current J R ). It is well known that an ax- 
ial current-carrying cylindrical plasma column with a helical 
magnetic field is subject to current-driven instabilities, such 
as sausage (m = 0), kink (rn = 1), and the other higher order 



m modes (m is the azimuthal mode number). We however do 
not see any visible evolution of the non-axisymmetric features 
in this magnetic tower jet. 

From this overview, we see that the magnetic tower jet can 
propagate through the stratified background medium while 
keeping well-ordered structures throughout the time evolu- 
tion. The magnetic fields push away the background gas, 
forming magnetically dominated, low-density cavities. This 
action also drives a HD shock wave which is ahead of and 
eventually separated from the magnetic structures. The mag- 
netic tower has a slender "body" from the confinement of the 
background pressure and an expanded "lobe" when the fields 
expand into a background with the decreasing pressure. 

We will now turn to the discussions on the detailed prop- 
erties of the tower jet, including the HD shock wav e and its 
imp act in the axial (z) and radial (x) directions in A3. 21 a nd 
13.31 The radial force balance of the jet is examined in A3 .41 

3.2. Structure of a Magnetic Tower Jet in the Axial (z) 
Direction 

Figure |3 displays several physical quantities along a line 
with (x, y) — (1, 0) in the axial direction at t = 0.75. Sev- 
eral features can be identified. First, the HD shock wave front 
can be seen around z ~ 13.5 in the profiles of p (top panel), 
V z , and C s (bottom panel). C s is higher than the initial back- 
ground value 1.29 in the post-shock region due to shock heat- 
ing and becomes smaller than 1.29 at z ~ 10.7 due to axial 
expansion. V z has the similar behavior as C s . 

Second, a magnetic tower front ("tower front" in the fol- 
lowing discussions) is located at z ~ 8.0, beyond which the 
magnetic field goes to zero, as seen in the top panel. This 
indicates that the gas within the magnetic tower jet is sepa- 
rated from the non-magnetized ambient gas beyond the tower 
front. We regard this front as a tangential discontinuity as the 
magnetic fields are tangential to the front without the normal 
component. This is consistent with the fact that the radial 
and azimuthal field components (B r and B^) are dominant 
near the tower front (z < 8.0) but the axial field component 
B z becomes dominant only for z < 6.0. The density and 
pressure show smooth transition through this front though the 
gradients of p and C s are slightly changed there. 

Third, there is another MHD wave front at z ~ 7.0 where 
B r , C s , and every velocity component have their local max- 
ima (p instead has its local minimum), as seen in both panels. 
To better understand the nature of this MHD wave front, we 
plot the axial profiles of pressures and various forces along 
the line (x, y) — (1, 0) at t — 7.5 in Fig. |6] The total pres- 
sure ptot consists of only the gas pressure p beyond the tower 
front (z > 8.0) but is dominated by the magnetic pressure 
p m behind the MHD wave front (z < 7.0), as seen in the top 
panel. A transition occurs around 7.0 < z < 8.0, where an 
increase in p is accompanied by a decrease in p m . We there- 
fore identify this as a reverse slow-mode compressional MHD 
wave front. In magnetic towers, the transition region between 
gas and magnetic pressures can be identified as a reverse slow 
wave front in the context of MHD wave structures. It does 
not depend on the resolution and parameters. So, the reverse 
slow mode wave (sometimes, it can be steepen into a shock) 
will always be there. In addition, in Fig. [7] we show several 
snapshots of V z and p during t = 7.5 ~ 10.0 (along the same 
offset axial path with Figs. |3]and[6j. The axial flow is deceler- 
ated by the gravitational force in the post-shock region beyond 
the tower front as seen in the top panel. On the other hand, 
the narrow region between the tower front and the reverse 
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FIG. 1. — The density evolution caused by a magnetic tower jet propagating through a stratified background. Color images and contours of the density 
distribution (logarithmic scale) are shown along with the poloidal velocity field (arrows) in the x — z plane. Density cavities are formed due to the magnetic field 
expansion. Times are given at the upper left in each of the four panels. The length of the arrow at the upper right in each panel shows the unit velocity. 

steepen into the reverse slow-mode MHD shock wave via this 
nonlinear evolution. Consequently, in the frame co-moving 
with the reverse slow-mode shock, a strong compression oc- 
curs behind the shock wave front and causes a local heating, 
as seen in the bottom panel and also Fig. |2] This heating could 
have interesting implications for the enhancement of radiation 
from radio to X-rays at the terminal part of Fanaroff-Riley 
type II AGN jets, such as lobes and hot spots, which are gen- 
erally interpreted as heati ng caused by the jet t erminal shock 
wave ( Blandfo rd~& Rees 1 1971 IScheuer 1 1 974 . 

Fourth, the HD shock wave breaks the initial background 
hydrostatic equilibrium. The passage of the shock wave heats 
the gas and alters its pressure gradient. As shown in the bot- 
tom panel of Fig. [6] the gas pressure gradient force F p stays 
uupositive at the shock front (which pushes the shock for- 
ward), but the total (gravity plus pressure gradient) force Ftot 
becomes negative behind the shock, implying a deceleration 
of the gas in the axial direction in the post-shock region. This 
is consistent with Fig. 
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FIG. 2. — The change of the gas pressure Ap/pi (Ap = p — p\) (pi 
represents p| t= o) at t = 10.0 in the x — z plane. 



slow-mode wave front is accelerated by the magnetic pres- 
sure gradient (the "magnetic piston" effect). Note that the re- 
verse slow-mode MHD compressional wave could eventually 



3.3. Deformation of the Jet "Body" in the Radial Direction 

We next examine the structure and dynamics of the mag- 
netic tower jet along the radial direction in the equatorial 
plane with (y, z) = (0, 0). Figure [8] shows the radial profiles 
of physical quantities along the x-axis at t = 6.0. The bound- 
ary of the magnetic tower jet ("tower edge" in the following 
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FIG. 3. — Three-dimensional configuration view of the selected magnetic 
lines at t = 1 0.0. 



discussions) is located at x ~ 3.0 where Fj x b becomes zero. 
Two distinct peaks of C s and V x around x ~ 7.8 and 9.5 are 
visible in the top panel. The first front (x ~ 9.5) is the prop- 
agating HD shock wave front as we showed in the previous 
section. The second front (x ~ 7.8) also indicates another 
expanding HD shock wave front generated by a bounce when 
the magnetic flux pinches in the radial direction caused by 
the "hoop stress". This secondary shock front appears only 
in the radial direction (see also Fig. 0. C s decreases gradu- 
ally towards the jet axis in the post-shock region of the sec- 
ondary shock and becomes smaller than its initial value 1.29 
at x ~ 5.7. We can confirm that these shock fronts are purely 
powered by the gas pressure gradients F p , as seen in the bot- 
tom panel (twin peaks of F p are shifted a bit behind that of C s 
in the top panel). 

Part of the region between the secondary HD shock front 
and the tower edge has F to t being negative (the bottom panel), 
meaning that the gas is undergoing gravitational contraction. 
This is indicated by V x < in the top panel for x < 5.1. This 



*= 10-0 




FIG. 4. — Distribution of the axial current density J z at t = 10.0 in the 
x — z plane. It shows a closed circulating current system, in which the current 
flows along the central (z) column ( J z > 0) and returns along the conically 
shaped shell on the outside ( J z < 0). 



behavior is similar to what we have discussed earlier, i.e., the 
HD shock waves break the background hydrostatic equilib- 
rium, causing a global contraction. Note that this contraction 
is occurring along the whole jet "body", e.g., for \z\ < 4.0 
when t — 7.5. 

The bottom panel also helps to address the question on 
what forces are confining the magnetic fields in the equato- 
rial plane. Since the total magnetic field Fj x b stays positive, 
this means that the inward hoop stress is not strong enough to 
confine the magnetic fields. Instead, at the tower edge, it is 
the background gravity that counters the combined effects of 
magnetic field pressure and a positive pressure gradient (push- 
ing outward). A bit further into the tower edge (at x ~ 2.1), 
however, the pressure force changes from outwardly directed 
to inwardly directed. Then, both the pressure gradient and 
the gravitational forces act to counter the outward J x B 
force. This behavior, which is mostly caused by the magnetic 
tower expanding in a background gravitational field, is differ- 
ent from the usual MHD models for jets where the inwardly 
directed hoop stress is bal anced by the outwardly d irected 
magnetic pressure gradient (Blandford & Pavne 1982). 

To investigate the radius evolution of the magnetic tower 
jet, we show the radial profile of density at the equatorial 
plane from t = 6.0 to 10.0 in Fig. [9] It shows that the 
radius has an approximately constant contraction speed for 
t = 6.5 ~ 9.5. The time scale for contraction is r contr ~ 6, 
which is about 7.5 times longer than the local sound-crossing 
time scale. 

3.4. Force Balance in the Radial Direction 

We now discuss the jet properties along the radial direc- 
tion away from the equatorial plane. Figure EH shows the 
radial profiles of physical quantities along the x-axis with 
(y, z) = (0, 4) at t = 7.5. The tower edge is now located 
at x ~ 3.0. The plasma j3 in the core of the tower is (3 < 0.1. 
The top panel shows that the central total pressure (which is 
dominated by the magnetic pressure) is much bigger than the 
background "confining" pr essure. This illustrates the original 
suggestion by Lynden-Bell (iLvnden-Bell 1 1 9961 12003 ) that the 
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FIG. 5. — Axial profiles of physical quantities, parallel to the 2-axis with 
(x, y) = (1, 0) at t = 7.5. Top: Density p and magnetic field compo- 
nents (B r , B^, B z ). Bottom: Sound speed C s and velocity components 
{V r , Vd, Vz). The positions of the expanding hydrodynamic shock wave 
front and the magnetic tower front, which is identified as a tangential dis- 
continuity, are shown in both panels. A horizontal solid line in the bottom 
panel denotes the initial sound speed (constant throughout the computational 
domain). 




FIG. 6. — Similar to Fig. |5] Top: Shown are the gas pressure p (dashed 
line), magnetic pressure (solid line) p m , and total pressure ptot (= P + Pm) 
(light gray thick solid line). Bottom: Shown are the forces in the axial (z) 

direction: the Lorentz force Fj x b 

line), gas pressure gradient force F p 



tational force F s 



JxB 



d/dz UBf +B^)/2j, (solid 

—dp/dz (dashed line), gravi- 
-pdtp/dR X \z\/R (dotted line), and total force 
- Fg) (light gray thick solid line). The position of the 



reverse slow-mode MHD wave front is also shown in the top panel. 



hoop stress of the toroidal field component B$ can act as a 
pressure amplifier in the central region of the magnetic tower: 
the pinch effect amplifies p m near the axis of the tower. At 
the tower edge, h owever, a fini te (albeit small) gas pressure is 
required (see also L i et al. 120011) . 

The bottom panel shows the detailed distributions of forces 
along the radial direction. They show that: 

• (Region: x > 3.0) Beyond the tower edge, the gravi- 
tational force F g is slightly stronger than the outwardly 
directed gas pressure gradient force F p , indicating that 
this edge is subject to the gravitational contraction, as 
discussed in the previous section. 

• (Region: 2.0 < x < 3.0) Interior to the tower edge, 
the Lorentz force Fj x b is dominated by the outwardly 
directed magnetic pressure gradient force F Pm , and it 
is also larger than the inwardly directed F p , indicat- 
ing that the outer shell of the magnetic shell should be 
expanding at the relatively higher z, in contrast to the 
equatorial region (z = 0) where the tower radius is con- 
tracting. The hoop stress Fh p plays a minor role in the 



force balance around this region. Note that J flows 
inside this region. 

• (Region: x < 2.0) Inside the jet "body", Contributions 
from both F Pm and Fh p become comparable and nearly 
cancel each other. The residual Fj x b is balanced by 
F p everywhere in this region. F^ p becomes dominant 
in the Lorentz force at the core part. Note that J F flows 
within x < 1.0. 

In addition, both the F g and the centrifugal force F c play 
a minor role in terms of the force balance inside the mag- 
netic tower jet. Thus, the interior of the magnetic tower 
is magnetically dominated but not exactly force-free, i.e., 
— Vp + J x B ~ . This small but finite pressure gra- 
dient force could potentially provide some stabilizing effects 
on the traditionally kink-unstable twisted magnetic configu- 
rations. The detailed examinations of stability properties in 
magnetic tower jets will be discussed in a forthcoming paper. 

4. DISCUSSIONS 
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FIG. 7. — Similar to Fig.|5] but with selected snapshots during t = 7.5 ~ 
10.0 (each time-interval is equal to 0.5). Top: The axial velocity V z . Bottom: 
The gas pressure p. 



On scales of ~ tens of kpc to hundreds of kpc, the back- 
ground density and pressure profiles have a strong influence 
on the overall morphology of the magnetic tower jet. Most no- 
tably, the transverse size of the magnetic tower grows as the 
jet propagates into a decreasing pressure environment, show- 
ing a similar morphology with the jet/lobe configuration of 
radio galaxies. 

The radial size of the lobe n be can be estimated from the 
following consideration: Figure fTTI shows the radial profiles 
of magnetic field components parallel to the x-axis with z = 
7.0 at t = 10.0. Together with Fig. g] we see that both the 
poloidal field and especially the poloidal current I z maintain 
well collimated around the central axis even at the late stage 
of the evolution. This implies that the toroidal magnetic fields 
in the lobe region are distributed roughly as fL~ I z /r. This 
is consistent with the results shown in Fig. II II where rB^ 
have a plateau from x m 1 — 2.5. As indicated in Fig. ^| we 
see that the magnetic pressure and the background pressure 
try to balance each other at the tower edge. Thus, we expect 
that 

Bl 

' 'Pe, (7) 



8tt 

where p c is the external gas pressure at the tower edge. 
When the lobe experiences sufficient expansion, we expect 
the poloidal field strength to drop much faster than l/r\ \, e . 





-d/dr[{Bl + Bl)/2] 



FIG. 8. — Radial profiles of physical quantities along the x-axis in the 
equatorial plane with (y, z) = (0, 0) at t = 6.0. Top: The sound speed 
G' B (light gray thick solid line) and the radial velocity V x (solid line). Bot- 
tom: The forces in the radial (x) direction: the Lorentz force Fj x g = 

B^/r (solid line), gas pressure gradient force 

'r (dashed line), gravitational force F g = — pdip / 8R X 
(dotted line), and total force Ftot (= FjxB + Fp + F g ) (light gray thick 
solid line). The positions of two expanding hydrodynamic shock wave fronts 
and the edge of the magnetic tower (tangential discontinuity) are shown in 
the top panel. A horizontal solid line in the top panel denotes the initial 
sound speed (constant throughout the computational domain) and a horizon- 
tal dashed line in the top panel represents a level "0" for the transverse veloc- 
ity. 



So we have 



which gives 



B % {Iz/riobc) 2 



8?r 



Stt 



n 



. - I T)- 1/2 
obc ± z Fe 



(8) 
(9) 



This is generally consistent with our numerical result shown 
in Fig. |l]thougn it is difficult to make a firm quantitative de- 
termination since the lobes have not expanded sufficiently. 

To make direct comparison between our simulations and 
observations, further analysis is clearly needed. Note that for 
the magnetic tower model, the Alfven surface is located at the 
outer edge of the magnetic tower and flow within the magnetic 
tower is always sub-Alfvenic. This is quite different from the 
hydromagnetic models where the MHD flow is accelerated 
and has passed through several critical velocity surfaces, in- 
cluding the Alfven surface. 
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FIG. 9. — The radial profiles of density p in the equatorial plane with 
selected snapshots during t = 6.0 ~ 10.0 (each time-interval is equal to 
0.5). The initial profile (t = 0.0) is also shown (dashed line). 

5. CONCLUSION 

By performing 3-D MHD simulations we have investigated 
in detail the nonlinear dynamics of magnetic tower jets, which 
propagate through the stably stratified background that is ini- 
tially in hydrostatic equil ibrium. Our s imulations, based on 
the approach developed in lLi et alTI J2006I) . confirm a number 
of the global characteristics developed in Lynden-Bell (1996, 
2003) and Li et al. (2001). The results presented here give 
additional details for a dynamically evolving jet in a strati- 
fied background. The magnetic tower is made of helical mag- 
netic fields, with poloidal flux and poloidal current concen- 
trated around the central axis. The "returning" portion of the 
poloidal flux and current is distributed on the outer shell of the 
magnetic tower. Together they form a self-contained system 
with magnetic fields and associated currents, being confined 
by the ambient pressure and gravity. The overall morphol- 
ogy exhibits a confined magnetic jet "body" with an expanded 
"lobe". The confinement of the "body" comes jointly from the 
external pressure and the gravity. The formation of the lobe is 
due to the expansion of magnetic fields in a decreasing back- 
ground pressure. 

A hydrodynamic shock wave initiated by the injection of 
magnetic energy/flux propagates ahead of the magnetic tower 
and can break the hydrostatic equilibrium of the ambient 
medium, causing a global gravitational contraction. As a re- 
sult, a strong compression occurs in the axial direction be- 
tween the magnetic tower front and the reverse slow-mode 
MHD shock wave front that follows. Furthermore, the mag- 
netic tower jets are deformed radially into a slender-shaped 
body due to the inward-directed flow of the ambient (non- 
magnetized) gas. 

The lobe is magnetically dominated and is likely filled with 
the toroidal magnetic fields generated by the central poloidal 
current. At the edge of the magnetic tower jet, the outward- 
directed magnetic pressure gradient force is balanced by the 
inward-directed gas pressure gradient force, so the radial 



width of the magnetic tower can be determined jointly by the 
magnitude of the poloidal current and the magnitude of the ex- 
ternal gas pressure. The highly wound helical magnetic field 
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FIG. 10. — The radial profiles of physical quantities along the x-axis with 
(y, z) = (0, 4) at t = 7.5. Top: Similar to the top panel in Fig. |5| Bot- 
tom: Similar to the bottom panel in Fig. |S] but with the centrifugal force 
F c = pV^ /r (dash-dotted line) and the Lorentz force Fj x g , which can be 

separated into the magnetic pressure gradient force —d/dr |A^^ + Bz )/^| 
(light gray thick solid line) and the hoop stress (magnetic tension force) 
—B^/r (dark gray thick solid line). The position of the magnetic tower 
edge (tangential discontinuity) is shown in both panels. 



in the magnetic tower never reaches the force-free equilibrium 
precisely, but obtains radial force-balance, including the gas 
pressure gradient inside the magnetic tower. 

The stability of the magnetic tower jets will be considered 
in our forthcoming papers. 
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TABLE 1 

Units of Physical Quantities for Normalization. 



Physical Quantities Description Normalization Units Typical Values 



R (= y/ x a + y 2 +z 2 ) Length R 5 Kpc 

V Velocity field C s0 4.6 X 10 7 cm/s 

t Time Ro/C s0 1.0 X 10 7 yrs 

p Density po 5.0 X 10 -27 g/cm 3 

p Pressure PoC 2 1.4 X 10 -11 dyn/cm 2 

B Magnetic field jAwp C 2 17.1 /iG 



NOTE. — The initial value of the density po and sound speed C s o at the origin [w, y, z) — (0, 0, 0) are chosen to be the typical density and velocity in the system. The initial 
dimensionless density p' and pressure p' at the origin are set to unity. We therefore have the initial dimensionless sound speed C' s — ~ 1.29. A characteristic time scale, the 

initial sound crossing time r s o — B-o / C s o — 10.0 Myrs, which corresponds to the dimensionless time scale r PS 0.78. 



t= 10.0 z= 7.0 
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FIG. 11. — The radial profiles of the magnetic field components 
(B r , B^, B z ), the poloidal magnetic field B p = J B 2 + B 2 , and the 
quantity rB^ along the rr-axis with (y, z) = (0, 7) at t = 10.0. The posi- 
tion of the magnetic tower edge (tangential discontinuity) is shown. 



Blandford, R. D., & Payne, D. G. 1982, MNRAS, 199, 883 

Churazov, E., Forman, W., Jones, C, & Bohringer, H. 2003, ApJ, 590, 225 

Ferrari, A. 1998, ARA&A, 36, 539 

Hsu, S. C, & Bellan, P. M. 2002, MNRAS, 334, 257 

Kato, Y., Hayashi, M. R., & Matsumoto, R. 2004, ApJ, 600, 338 

Kato, Y., Mineshige, S., & Shibata, K. 2004, ApJ, 605, 307 

King, I. 1962, AJ, 67, 471 

Kudoh, T., Matsumoto, R., & Shibata, K. 2002, PASJ, 54, 26 
Lebedev, S. V. et al. 2005, MNRAS, 361, 97 

Li, H„ Lovelace, R. V. E., Finn, J. M., & Colgate, S. A. 2001, ApJ, 561, 915 
Li, Ft., Lapent a, G„ F inn, J. M., Li, S., & Colgate, S. A. 2006, ApJ in press 
1 astro-ph/0604469 1 

Li, y., & Li, H. 2U03, Technical Report, LA-UR-03-8935, Los Alamos 

National Laboratory 
Lovelace, R. V. E. 1976, Nature, 262, 649 

Lovelace, R. V. E., Li, Ft., Koldoba, A. V., Ustyugova, G. V., & Romanova, 

M. M. 2002, ApJ, 572, 445 
Lovelace, R. V. E., & Romanova, M. M. 2003, ApJ, 596, L159 
Lynden-Bell, D. 1996, MNRAS, 279, 389 

Lynden-Bell, D. 2003, MNRAS, 341, 1360 

Lynden-Bell, D. 2006, MNRAS in press 1 astro-ph/0604424 1 

Lynden-Bell, D., & Boily, C. 1994, MNRAS, 267, 146 

Meier, D. L., Koide, S., & Uchida, Y. 2001, Science, 291, 84 

Romanova, M. M., Ustyugova, G. V., Koldoba, A. V., Chechetkin, V. M., & 

Lovelace, R. V. E. 1998, ApJ, 500, 703 
Scheuer, P. A. G. 1974, MNRAS, 166, 513 

Turner, N. J., Bodenheimer, P., & Rozyczka M. 1999, ApJ, 524, 12 
Uchida, Y, & Shibata, K. 1985, PASJ, 37, 515 

Ustyugova, G. V., Lovelace, R. V. E., Romanova, M. M., Li, H., Colgate, S. 

A. (2000), ApJ, 541, L21 
Uzdensky, D. A., & MacFadyen, A. I. 2006, preprint I astro-ph/0602419 1 



Blandford, R. D., & Znajek, R. L. 1977, MNRAS, 179, 433 



